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An extended Landau-Levich model of liquid-propellant combustion, one that allows for a local depen- 
dence of the burning rate on the (gas) pressure at the liquid-gas interface, exhibits not only the classical 
hydrodynamic cellular instability attributed to Landau but also a pulsating hydrodynamic instability asso- 
ciated with sufficiently negative pressure sensitivities. Exploiting the realistic limit of small values of the 
gas-to-liquid density ratio p , analytical formulas for both neutral stability boundaries may be obtained by 
expanding all quantities in appropriate powers of p in each of three distinguished wave-number regimes. 
In particular, composite analytical expressions are derived for the neutral stability boundaries A p (k), where 
A J? is the pressure sensitivity of the burning rate and k is the wave number of the disturbance. For the 
cellular boundary, the results demonstrate explicitly the stabilizing effect of gravity on long-wave distur- 
bances, the stabilizing effect of viscosity (both liquid and gas) and surface tension on short-wave pertur- 
bations, and the instability associated with intermediate wave numbers for negative values of A p> which is 
characteristic of many hydroxyiarnmonium nitrate-based liquid propellants over certain pressure ranges. 
In contrast, the pulsating hydrodynamic stability boundary is insensitive to gravitational and surface-tension 
effects but is more sensitive to the effects of liquid viscosity because, for typical nonzero values of the 
latter, the pulsating boundary decreases to larger negative values of A p as k increases through 0(1) values. 
Thus, liquid-propellant combustion is predicted to be stable (that is, steady and planar) only for a range 
of negative pressure sensitivities that he below the cellular boundary that exists for sufficiently small 
negative values of A p and above the pulsating boundary that exists for larger negative values of this param- 
eter. 


Introduction 

The stability of liquid-propellant deflagration is a 
fundamental problem that was first treated by Lan- 
dau in a classical study [1] that introduced the con- 
cept of hydrodynamic instability in a combustion 
context. Referred to as the “slow combustion of liq- 
uids/’ that analysis is most applicable to certain re- 
alistic, limiting cases in which combustion may be 
approximated by an overall reaction at the liquid- 
gas interface. For example, the gas flame may occur 
under near-breakaway conditions, exerting little 
thermal or hydrodynamic influence on the burning 
propellant; or distributed combustion may only oc- 
cur in an intrusive regime such that the reaction 
zone lies closer to the liquid-gas interface than the 
length scale of any disturbance of interest; or the 
liquid propellant may simply undergo exothermic 
decomposition at the surface without any significant 
distributed combustion, such as appears to occur in 
some types of hydroxyiarnmonium nitrate (HAN)- 
based liquid propellants at low pressures [2]. The 


results of Landau’s study, along with a subsequent 
paper by Levich [3] that replaced the effects due to 
surface tension in the earlier study with those due 
to (liquid) viscosity, have been widely quoted and 
offer much in terms of physical insight into the na- 
ture of this type of instability, which, as in the case 
of gaseous combustion, is associated with the density 
change across the reaction front. However, because 
these models assumed a constant normal burning 
rate, it has proven useful to improve upon them by 
incorporating a more realistic coupling of the normal 
propagation speed with the local pressure and tem- 
perature fields, thereby allowing for a locally varying 
burning rate [4,5]. One result that has emerged from 
this generalization is that, in addition to the classical 
Landau (cellular) type of instability, the models now 
predict a pulsating hydrodynamic instability as well. 
The latter arises specifically from the local pressure 
coupling and thus may be physically achievable be- 
cause the mass burning rate of many propellants has 
been shown empirically to correlate well with 
pressure. Thermal coupling, on the other hand, 
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introduces additional thermal/diffusive instabilities 
[5] that will not be considered here. In the present 
work, we shall consider both types of hydrodynamic 
instabilities, but particular focus will be placed on 
the pulsating stability boundary, because this type of 
instability is absent from the earlier models that ne- 
glected the pressure coupling previously indicated, 
As in a companion study that focused on the classical 
(cellular) hydrodynamic instability [6], we develop a 
formal asymptotic theory by considering the realistic 
limiting case in which the gas-to-liquid density ratio 
p is small. 
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The governing hydrodynamic equations consist of 
mass and momentum on Either side of the gas-liquid 
interface, supplemented by a pressure-dependent 
burning-rate law and associated continuity and jump 
conditions across the interface [4—6]. Thus, it is as- 
sumed, as in the classical models, that there is no 
distributed reaction in either the liquid or gas phases 
but that there exists either a pyrolysis reaction or an 
exothermic decomposition at the liquid-gas inter- 
face that depends on the local pressure. For sim- 
plicity, it is assumed that within the liquid and gas 
phases separately, the density and other fluid prop- 
erties are constants, with appropriate jumps across 
the phase boundary. The nondimensional location of 
the fatter is denoted by x 3 = & s (x lt x 2 , *)> where the 
adopted coordinate system is fixed with respect to 
the stationary liquid at x 3 = — °°. Then, in the mov- 
ing coordinate system x = x l5 y = x 2 , z = x 3 - 
<p s (x h x 2 , t), in terms of which the liquid-gas inter- 
face always lies at z = 0, the complete nondimen- 
sional formulation of the problem in the absence of 
thermal coupling is given by 

V v = 0, z ^ 0 (1) 

dv _ dv + (v . v)v = (0 o _ Fr -i) 

at at dz 
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ditions 
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where the latter represent continuity of transverse 
velocity components (no slip), conservation of (nor- 
mal) mass flux, the mass burning-rate law, and con- 
servation of normal and transverse components of 
momentum flux, respectively. Here, v and p denote 
velocity [with respect to the original (x t , x 2 , x 3 ) co- 
ordinate system] and pressure, and the ± subscripts 
denote evaluation at z = 0*. The parameters Pr, and 
Pr g denote the liquid and gas-phase Prandtl num- 
bers, p and l are the gas-to-liquid density and ther- 
mal diffusivity ratios, Fr is the Froude number, e is 
the rate-of-strain tensor, and y is the surface-tension 
coefficient. All of these quantities are defined in 
terms of their dimensional counterparts in the no- 
menclature. In addition, the factor S(^. v ) and the 
unit normal ri s are given by = [1 + (B<P s /Bx) 2 
+ (d<P s /dy) 2 ]- l/2 and n s = (-30,/ax, -B& s /By y 
1 )S(<£ V ), while the expressions for the gradient 
operator V and the Laplacian V 2 in the moving 
coordinate system are given by V = [d/dx — 

(, d0 s !dx)d/dz , Sidy ~ (B& s /dy)d/Bz, d/dz] and 
V 2 = B 2 /Bx 2 4 &/By 2 + [1 + (B&Jdx) 2 4 
(d<P s /< dy) 2 ]B 2 /dz 2 - 2 (d& s /dx)d 2 /dxdz 

2(d<P s /dy)B 2 /dyBz - (d 2 0 s /dx 2 4 d 2 & s /By 2 )d/dz. We 
remark that the factor multiplying y in equation 4 is 
the curvature — V ■ A, and note that pAPr g = 
where p is the gas-to-liquid viscosity ratio. 

Finally we observe that the burning rate A(p) in 
the last of equations 3 is assumed, in accordance with 
numerous experimental correlations for both solid 
and liquid propellants [7], to depend on the local gas 
pressure at the interface, where A is normalized to 
unity for the case of steady, planar burning. Indeed, 
in the linear stability analysis that follows, the pres- 
sure sensitivity A p = BA/dp\ p „ 0 (where, as indicated 
in equation 6 in the following, p — 0 is the unper- 
turbed gas pressure at z = 0 + ) of the local burning 
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rate will emerge as an important parameter. In fact, 
A p is the only information about A (p ) that will appear 
explicitly, although in a nonlinear stability analysis, 
more detailed information in the form of higher de- 
rivatives would enter into the calculations (cf. Ref. 
[8]). We remark that, depending on the particular 
oropellant, A p is not necessarily positive and may in 
fact take on zero and negative values over certain 
pressure ranges [2,7]. 

A nontrivial basic solution to the foregoing prob- 
lem, corresponding to the special case of a steady, 
planar deflagration, is given by 

<&> = -t, v° = (0, 0, v°) 




C(x, ij , z , t) are defined as + <f> s > v = 

v°(z) + u, and p = p°(z) + C, respectively. Substi- 
tuting these definitions into the nonlinear model de- 
fined earlier and linearizing about the basic solution, 
the perturbation problem becomes 
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The linear stability analysis of this solution now pro- 
ceeds in a standard fashion. However, owing to the 
significant number of parameters, a complete anal- 
ysis of the resulting dispersion relation is quite com- 
plex, and we follow our previous approach [6] by 
restricting further consideration to the realistic pa- 
rameter regime 1,|U < 1 and, in the case of 
microgravity, Fr~ 1 C 1. In contrast, the earlier clas- 
sical studies only considered special limiting cases 
and/or assumptions. Thus, in the study due to Lan- 
dau [1], viscosity was neglected and the effects of 
gravity (assumed to act normal to the undisturbed 
planar interface in the direction of the unbumed liq- 
uid) and surface tension were shown to be stabiliz- 
ing, leading to a criterion for the absolute stability 
for steady, planar deflagration of the form (in our 
nondimensional notation) AyFr~ x fP/(\ — p) > 1. In 
the study due to Levich [3], surface tension was ne- 
glected, but the effects due to the viscosity of the 
liquid were included, leading to the absolute stability 
criterion Fr~ 1 Fr / ( 3/7 ) V2 > 1. Thus, these two studies, 
both of which assumed a constant normal burning 
rate (A = 1), demonstrated that sufficiently large 
values of either viscosity or surface tension, when 
coupled with the effects due to gravity, may render 
steady, planar deflagration stable to hydrodynamic 
disturbances. In our recent study [6], these results 
were synthesized and extended to the more realistic 
case of a nonconstant burning rate (i.e., A p ^ 0) in 
the limiting parameter regime identified earlier. In 
the present work, we summarize these results for the 
classical cellular boundary and use the resulting scal- 
ings to derive an expression for the pulsating hydro- 
dynamic stability boundary that arises from the pres- 
sure dependence of the local burning rate. 

Linear Stability Problem 

With respect to the basic solution, equation 6, the 
perturbation quantities <k(%, y , 0, u(x, y, z, t), and 
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equations 11-13. 

Nontrivial harmonic solutions for u, and £, pro- 
portional to ef °* +ik ' that satisfy equations 1 and 
2 and the boundary and boundedness conditions at 
z — ± co axe given by 

^ _ e i(o( +ik ] x + ik i y 


f b.e^ - Fr-\ z < 0 

< = +*•**+** \ b ^. kz . pFr . K z > o 
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where the foregoing solution has been normalized 
by setting the coefficient of the harmonic depen- 
dence of 4>, to unity. Here, the signs of k , and k, 
may be either positive or negative, and we have em- 
ployed the definitions k = (k ? + k\) m , and q and r 
are defined as 2 Prfl = 1 + [1 + 4 Pr,(iw + 
Pr,k 2 )] v2 and IpPrf = 1 - [1 + 4pPr,(imp + 
pPr t k 2 )] m . 

Substituting this solution into the interlace con- 
ditions 9-13 and using equation 7 for z ^ 0 yields 
nine conditions for the eight coefficients bi-b$ and 
the complex frequency (dispersion relation) ico(k). 
In particular, these conditions are given by 
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Although the foregoing problem is linear in the co- 
efficients b x -b H , explicit expressions for the disper- 
sion relation ico(k) and the neutral stability bound- 
aries are not readily obtainable in closed form, 
except in certain special cases as noted later. How- 
ever, it is possible to develop tractable perturbation 
expansions for these quantities in the realistic limit 
that the density and viscosity ratios p and// are small, 
as is Fr" 1 in the case of reduced gravity. 


Asymptotic Analysis of the Cellular Stability 
Boundary 

It turns out that the p C 1 limit implies the exis- 
tence of several different wave-number regimes [6], 
which in turn implies several different expansions for 
the dispersion relation. This is motivated by consid- 
ering the solution of equations 18-25 in the limit of 
zero viscosity (// = = 0), which leads to a trac- 

table form of the dispersion relation for arbitrary /?. 
In particular, the neutral stability boundaries with 
respect to infinitesimal hydrodynamic disturbances 
proportional to e l “** ik ' x , where k and x are the trans- 
verse wave number and coordinate vectors, respec- 
tively, are given for A p < 0 by [5] 


p( 1 - p)Fr ~ ] 4- pyk 2 - (1 “ p)k 
P p \ 3 - p)Fr~ x + p 2 yk 2 + (1 - />)(2 - p)k 

< 0 , co = 0 

(26) 

and 

A p = -p/( 1 - P), a? = k + P ) 

.Q _ pj-ipr" 1 + k 2 ( 1 - p)~ l y + kJp] (27) 

where k = lid. For A p > 0, the basic solution is 
always unstable, and thus steady, planar combustion 
is only stable in the region A p < 0 that lies between 
these two curves (Fig. 1). The first of these bound- 
aries is a cellular boundary (or = 0) that corresponds 
to the classical Landau instability. Indeed, in the 
limit A = 0, corresponding to a constant normal 
burning rate independent of pressure, Landau s clas- 
sical result is recovered since in the limit that yFr ~ 
approaches the value (1 - p)/4p 2 from below, the 
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Fig. 1. Hydrodynamic neutral sta- 
bility boundaries in the limit of zero 
viscosity. 


cellular stability boundary recedes from the region 
A p < 0. The pulsating stability boundary (27), on the 
other hand, only occurs for nonzero Ap and, hence, 
was not predicted by the classical theories that as- 
sumed a constant normal burning rate. We note that 
zero and negative values of A p over certain pressure 
ranges are characteristic of the HAN -based liquid 
propellants mentioned earlier [2]. 

The formalism necessary to analyze the fully vis- 
cous problem is suggested by the fact that, for small 
p 7 the cellular boundary (26) has different limiting 
forms dependent on the relative magnitude of A: with 
respect to p . Thus, based on characteristic parameter 
values, we introduce a bookkeeping parameter e C 
1 and define the scaled parameters g*, p*, p*, and 
A* according to p — p *£ , p = //*£, A p = A*£, and 
either Fr~ l = g or Fr~ l = g*e, where y and Pr } s 
P are regarded as 0(1) and the scaling for A p is mo- 
tivated by equation 26. Here, the second scaling for 


Fr~ l corresponds to a reduced gravity limit, whereas 
the first definition indicates the normal gravity case. 
Equation 26 then suggests three wave-number 
scales; an inner (small) scale A,-, the outer 0(1) scale 
A, and a far outer (large) scale kr = Ae, where the 
inner scale is defined as A* = kfe for Fr 1 ~ 0(1) 
and k t = A/e 2 for the reduced gravity limit Fr' 1 ~ 
0(e). Thus, in each of these regions, the correspond- 
ing leading-order expressions for A^ are deduced 
from equation 26 as 




p*(p+ g - Km 

p*(p*g* - *i)/2A| 


A*<«) I P*, tfP ~ I p*(p*ykf - 1) (28) 


and hence a uniformly valid composite expansion 
A£( c) (k) may be constructed as 
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where the definitions of ^ &nd kj~ have been used to 
express the final result in terms of k , and the upper 
and lower expressions for A*( f) correspond respec- 
tively to the normal and reduced gravity limits de- 
scribed earlier. In terms of the original unsealed pa- 
rameters, equation 29 becomes 2^ p + yp 2 k 

+ Yr~ 1 p 2 /k, which, in the parameter regime consid- 
ered, is a leading-order asymptotic representation of 
the exact relation 26. It is readily deduced from 
equations 28 and 29, as discussed in further detail 
elsewhere [6], that surface tension stabilizes large 
wave-number disturbances, while gravity stabilizes 
small wave-number perturbations. In the reduced- 
gravity limit, the minimum in the stability boundary 
is thus shifted to smaller wave numbers, and thus 
the hydrodynamic cellular instability becomes a 
long-wave instability phenomenon (Fig. 2). 
Corresponding results may be obtained for the vis 
cous case. However, to deduce the asymptotic 
forms of both the cellular and the pulsating stability 
boundaries, it is preferable from the standpoint of 
tractability to introduce the scalings introduced pre- 
viously and appropriate perturbation expansions for 
the coefficients b t directly into equations 18-25 for 
each wave-number regime, and to obtain the neutral 
stability boundary in question from the expanded 
form of those equations. Thus, for the cellular 
boundary A*(k) [6], we deduce the following expan- 
sions and results for each wave-number regime. 
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(i) k ~ 0(1) 

r ~ r x e + • • • , r, = p*Pk 2 
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q 0 = (2 P)- 1 [1 + (1 + 4 P 2 fc2)i«] (30) 
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We observe that the leading-order results 32 and 
35 are equivalent to the corresponding inviscid re- 
sults 28. Thus, to leading order, neither the inner 
nor the outer wave-number regimes are influenced 
by viscous effects, which, to a first approximation, 
are only significant for large wave-number distur- 
bances. This is reflected in the leading-order ex- 
pression for the cellular stability boundary given by 
equation 38, where, among other features, it is read- 
ily observed that both the liquid and the gas-phase 
viscosities (through the parameters P and p*P, re- 
spectively) enter into this expression, reflecting an 
equal influence of viscous and surface-tension ef- 
fects on the neutral stability boundary in the large 
wave-number regime. The equal importance of gas- 
phase viscosity relative to that of the liquid phase 
stems from the fact that gas-phase disturbances are, 
according to equations 37, larger in magnitude than 
those in the liquid phase, such that a weak damping 
of a larger magnitude disturbance is as significant as 
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HYDRODYNAMIC STABILITY BOUNDARIES (p « 1) 



Fig. 2. Asymptotic representation 
of the cellular hydrodynamic neutral 
stability boundary in the limit of zero 
viscosity. The upper {lower) solid 
curves correspond to the two cases 
described by equations 29 for normal 
and reduced-gravity regimes, re- 
spectively (curves drawn for the 
case £ = 0.04, p* = 1-0, g = 2.5, 

g* = i-o). 


an 0(1) damping of a smaller magnitude distur- 
bance. In the limit P -> 0, the inviscid expression 28 
is recovered. It is easily shown that lim^ “ 
- p*/ 2, so that the far outer solution can oe matched 
to the outer solution 32. Indeed, a uniformly valid 
composite expansion spanning all three wave-num- 
ber regimes may be constructed as in the inviscid 
case, giving the result 

AJW - -p* + 

2p*p*P[\ + ek(p*y + 2 p*P + 2p*P) ) 

4p*P( 1 + ekp'P) ~ ( p*y + - (1 + 4p*WW*} 


(39) 


Ajf > p + 

lppP[\ + k{py + tyP + 2pf)3 

4pP{l + kpP) - (py + 2pP)[l - (1 + 1 


(40) 

The cellular stability boundaries, based on equation 
39, are reproduced in Fig. 3, where only those por- 
tions of the curves that lie in the region A * ^ 0 are 
shown. For sufficiently small positive values of A p , it 
may be shown that there always exists a positive 
(real) root ia) of the dispersion relation, which im- 
plies that this region is intrinsically unstable. We 
note from equation 38 that as fcy increases, Aq ^ ^ in- 
creases, intersecting the A* (A — 0 axis at the value 
k f - (p*y)~ l [l - p*P/(p*y + 2p*P)], which agrees 


or, reverting to unsealed parameters. 
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Fig. 3. Asymptotic representation 
of the cellular hydrodynamic neutral 
stability boundary for the viscous 
case. The upper and lower sets of 
curves correspond to the normal and 
reduced-gravity regimes, respec- 
tively, in the asymptotic limit consid- 
ered in this work (curves drawn for 
the case £ = = 1.0, g = 6.0, 

g* = 2.0). The solid curves corre- 
spond to the inviscid limit (P = 0) 
with nonzero surface tension (y = 
2.5). The dash-dot curves corre- 
spond to nonzero surface tension (y 
= 2.5) and liquid viscosity (P = 1.0) 
but zero gas-phase viscosity ~ 
0). The dash-dot-dot curves differ 
from the dash-dot curves by the ad- 
dition of gas-phase viscosity (p*P = 
1.0) and are similar to the dash-dot- 
dot-dot curves, where the latter cor- 
respond to larger viscosities (P = 
- 2.0). The dash-dot-dot-dot- 
dot curves correspond to a viscous 
case (P = fi*P = 10) but with zero 
surface tension (y = 0), so that, from 
equation 38, the curves do not inter- 
cept the A* = 0 axis. 


with the inviscid result in the limit fi*P — » 0. It is 
readily seen from Fig. 3 that, as in the inviscid case, 
the essential qualitative difference between the 
normal and reduced-gravity curves is the location of 
the critical wave number for instability. Specifically, 
the minimum in the neutral stability boundaries oc- 
curs for 0(1) values of k under normal gravity, and 
at k — 0(e 1/2 ) in the reduced-gravity limit considered 
here. It is also clear from Fig. 3 that increasing the 
values of any of the parameters P } //*P, or y serves 
to shrink the size of the unstable domain through 
damping of short-wave perturbations. The non-neg- 
ligible effects of gas-phase viscosity represents an 
important correction to Levichs original treatment 
[3] in which these effects were simply assumed to 
be small. The results 39-40 thus synthesize and sig- 
nificantly extend the classical Landau-Levich results 
[1,3], not only in allowing for a dynamic dependence 
of the burning rate on local conditions in the vicinity 


of the liquid-gas interface but also in its formal treat- 
ment of those processes (surface tension, liquid and 
gas-phase viscosity) that affect damping of large 
wave-number disturbances. 

Asymptotic Analysis of the Pulsating Stability 
Boundary 

As indicated previously, the existence of a nonsta- 
tionary pressure dependence on the burning rate 
(i.e., Ap 7* 0 leads to the prediction of a pulsating 
hydrodynamic stability boundary that is absent when 
such a pressure coupling is neglected, as in the origi- 
nal Landau-Levich theories. In the inviscid case, 
this boundary (equation 27) is a straight line that lies 
below the cellular boundary discussed earlier, but 
this is modified under the influence of viscosity, as 
we shall demonstrate. 

For the scalings adopted in the preceding section 
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[in particular, for P ~ 0(1), p ~ O(c)], it turns out 
that, unlike the cellular stability boundary for which 
viscous effects only have a leading-order effect in the 
far outer wave-number regime, the effects of viscos- 
ity have a leading-order effect on the pulsating 
boundary for 0(1) wave numbers as well. Thus, in 
the outer wave-number region, we seek a solution 
for the dispersion relation in the form ico ~ s 1/2 (i&>o 
+ i(o } £ m + ia) 2 e l/2 + •**), where the leading-order 
tenn is suggested by the explicit results for the in- 
viscid case [5,6], and the expansion in powers off 
is suggested by the expansions for r and q given be- 
low equation 17, whicn have the form 

r - r (!/2) £ 1/2 + r m ) e y4 + r,« + • ■ ■ 
r<i/2) - -ico t) p*, r ( 3 / 4) = 
r, = -to** - (**»>* ( 41 > 

q ~ 9(_ 1/4)6 1/4 + < 7 o£° + 

<7, - 1/4) = (taVP)“ 

q Q = (2 P)' 1 [1 + iw l /(ico 0 /P) m ] (42) 

Corresponding expansions for the coefficients b. in 
equations 18-25 are determined as 

b , = + b\-™' 

+ £(-1/2) £-1/2 + i = 1, 2, 8 (43) 

i 

b, = £ 1/2 + £ ~ 1/4 + ■ 

i = 3, 4, 5, 6; 

b t = £,(-m) f -L« + fc| 0) £° +■■•,< = V (44) 

where the leading terms in the expansions for b\, hi, 
b 4 , b 6 , and b H are consistent with the inviscid results 
[5] and the remaining coefficients appear only for 
nonzero values of P and are conservatively postu- 
lated to have the indicated expansions. Substituting 
these expansions into equations 18—25 and equating 
coefficients of like powers of e, we obtain the lead- 
ing-order equations and results 

ik, + ik 2 bif'V + q ( _ 1/4) b*T 1/4> - 0 

= ik[ b [ r m) + ik 2 b\r m) + r (M > *8 1 

(45) 

b[~» + bif" - 2b^-J> = 0 
h(,-D = -kip*, b [~ 1 1 = (to 0 ) 2 /fc 

fo<-» = -(k/p*)(l + App*) (46) 

where the last of equations 46 was obtained from 
the leading-order difference of equations 21 and 22, 
and the remainder of the leading-order versions of 


equations 18-25 give redundant results. Combining 
equations 46, we obtain 

(to,,) 2 = (kVp*)( 1 + 2 A*/p*) (47) 

and thus (to,,) 2 « 0 for A* 5 -pV % which essen- 
tially recovers the leading-order cellular stability 
boundary (32) for 0(1) wave numbers, but gives no 
information on the pulsating boundary because to o 
is purely imaginary for A* < — p*/2. Hence, stabil- 
ity in the latter region is determined by higher-order 
coefficients in the foregoing expansion for to. 

Continuing with the analysis of the expanded 
forms of equations 18—25, we obtain the second-or- 
der equations and results 

il t, + ik 2 + q ( -m hf 

+ qc.H ' ,/4)=0 (48) 

it] b'f ,/4 ' + ikjy'f 1/4> + r {m) b\f M) 

+ r ( 3 / 4 ) bp l) = 0 (49) 

k - <?(- i/ 4 M' ,,4) = iWc (1 - A;/p*) 

bif 1/41 - (k/io>o)b[- m = to, ( 5 °) 

£,(-3/4) _ = b ( H m) = bif 1/2) 

= 1/2) = 0 (51) 

where the first of equations 50 was obtained from 
the sum of equation 19 multiplied by ik , and equa- 
tion 20 multiplied by ik 2 and the first of equations 
45 , the fact that b ( f m) = 0 follows from the differ- 
ence of equations 21 and 22 , and the last two of 
equations 51 follow from equations 24 and 25 in con- 
junction with the result bif^ = 0. From these re- 
sults and the first of equations 45, we thus conclude 
that 

bir i /4 > = t<»i = 0, 

b'f = to 0 (1 - App*) (52) 

where the fact that to, = 0 implies the need to 
calculate to 2 to determine stability in the region 
A* < - p * 72 . Proceeding in this fashion, we obtain 
from the next-order versions of equations 21 and 23, 
the difference of equations 21 and 22 , and the sum 
of equation 24 multiplied by ik, and equation 25 
multiplied by ik 2 , the relations 

bL 0) — (k/ia>fb { f im + 

[<to* + kWco n m b[-» 

- p*bf 1 1 + P*bif 1) = to 2 (53) 

b \- 1/2 > + (k/m,)(2Pk - l)bp l) + bif 1/2) 

+ C2ia>nP* !k)bif l) - 2b<f m) = -to 0 (54) 
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Fig. 4. Asymptotic representation 
of the pulsating hydrodynamic neu- 
tral stability boundary for the viscous 
case (P > 0). The region between the 
pulsating and cellular boundaries 
(the latter are shown on an expanded 
scale in Fig. 3) is the stable region 
with respect to hydrodynamic insta- 
bility. 


-p*b { H xm + p*¥ 2 ~ m + (im 0 p* 2 lk)bi- l) 

+ A*U 2 xm = - ico oP * ( 55 ) 

+ k 1/21 + i(o nP *b‘f' ) 

+ too*,!®’ - (2Pk - 1)(J WkOo)b\- l) = 0 (56) 

which when combined with the expressions for 
b^'\ b< H l \ tt 2 m \ and icon given earlier, con- 
stitute four equations for the four unknowns 
b[~ l/2 \ b^\ b^~ l/2 \ and to> 2 - Solving these simulta- 
neous equations, we thus obtain 

M?> = -2 Pk 2 , 

b (-i/2) _ [2pk + 1 + A */p* - 2(A*/p*) 2 )i<o 0 
b^ l/2) = [1 - (A*//)*) 2 ]tcuo 


ico 2 = k{(App*f - 2 Pk - 1] (58) 

Equation 58 is the desired result, from which we 
conclude that ia >2 5 0 f° r (A*/p* ) 2 $ 1 + 2 Pk. 
Thus, in the region A* < 0, ia> 2 vanishes on the 
boundary 


A* p*( 1 + 2Pfc) 1/2 (59) 

which is a pulsating boundary (Fig. 4) because, from 
equation 47, ia) 0 is purely imaginary along this curve. 

Equation 59 is valid for 0(1) wave numbers, but 
because it matches to the leading-order inviscid in- 
ner pulsating boundary A* ~ — p* as k — » 0 and 
becomes large in a negative sense as k becomes 
large, it is clear that equation 59 represents the pul- 
sating boundary for arbitraiy wave numbers. That is, 
for P ~ 0(1), the effects ot (liquid) viscosity on the 
pulsating boundary are, to a first approximation, ab- 
sent for small wave numbers, are first felt for 0(1) 
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wave-number perturbations, and are sufficient to 
move this boundary to larger-magnitude values 
A* ~ 0(e' 1/2 ) in the far outer wave-number re- 
gime. In contrast, the cellular boundary 39 is unaf- 
fected for 0(1) and smaller wave numbers and is 
only modified an 0(1) amount for 0(e~ l ) wave num- 
bers. Thus, the hydrodynamic pulsating boundary is 
more sensitive to viscous effects than is the corre- 
sponding cellular stability boundary. For smaller- 
magnitude viscosities such that P = Pe ~ O(e), it 
may be shown by an analogous calculation (Appen- 
dix) that 0(1) modifications to the pulsating bound- 
ary then occur in the far outer wave-number regime 
according to A* ~ — p*(l + 2 which, in 

terms of unsealed quantities, is the same as equation 
59. 


Conclusion 

The present work has presented a formal asymp- 
totic treatment of hydrodynamic instability for a sur- 
face model of liquid-propellant combustion in which 
burning takes place at the liquid-gas interface. The 
model itself is based on a synthesized version of the 
classical models analyzed by Landau [1] and Levich 
[3], generalized to allow a coupling of the burning 
rate with the local pressure field [4,5]. The realistic 
smallness of the gas-to-liquid density ratio proved to 
be a convenient small parameter upon which to base 
an asymptotic treatment, resulting in three distinct 
wave-number regimes with different physical pro- 
cesses assuming dominance in each. Both cellular 
and pulsating hydrodynamic stability boundaries are 
predicted by the present model, the former corre- 
sponding to Landaus original notion of hydrody- 
namic instability and the latter representing a new 
prediction arising from the pressure dependence of 
the burning rate. For the cellular type of instability, 
it was shown that the gravitational acceleration (as- 
sumed to be normal to the undisturbed liquid-gas 
interface in the direction of the liquid) is responsible 
for stabilizing long- wave disturbances, whereas sur- 
face tension and viscosity are effective in stabilizing 
short-wave perturbations. In the case of pulsating 
instability, neither gravity nor surface tension play a 
leading-order role, and viscous effects are the dom- 
inant stabilizing influence. Indeed, for 0(1) liquid 
Prandtl numbers, the stabilizing effects of (liquid) 
viscosity on pulsating instability are significant for 
disturbances whose wave numbers are 0(1) and 
higher. On the other hand, viscous effects are only 
significant for large wave-number disturbances in 
the case of cellular instability, where the influence 
of gas and liquid viscosity are comparable despite the 
small ratio of these two parameters. Although the 
onset of pulsating hydrodynamic instability is pre- 
dicted to occur only for sufficiently negative values 


of the pressure-sensitivity coefficient A }} , the persis- 
tence of the pulsating stability boundary (in the pres- 
ence of viscous effects) for small wave numbers sug- 
gests that it should be observable in those types of 
liquid propellants, such as those based on hydroxy- 
lammonium nitrate (HAN) and triethanolammon- 
ium nitrate (TEAN), that are characterized by neg- 
ative pressure sensitivities over certain pressure 
ranges. In connection with this, we note that sloshing 
behavior has been observed during combustion of 
certain HAN -TEAN -water mixtures [2], but because 
nonsteady burning can arise via secondary and 
higher-order bifurcations in the cellular region [8], 
as well as from a primary crossing of the pulsating 
boundary described here, further measurements are 
generally needed to determine the precise origin of 
such behavior in any given experiment. 


Nomenclature 

A burning rate 

A p pressure-sensitivity coefficient 

b i coefficients in perturbation solution ( i = 

1, 2, . . . , 8 ) 

e rate- of- strain tensor 

Fr Froude number 

g inverse Froude number (gravitational ac- 

celeration) 

k perturbation wave number 

ft s unit normal 

p pressure 

P t Pr Prandtl number 

q quantity defined below equation ( 17) 

r quantity defined below equation ( 17) 

t time variable 

u perturbation velocity vector 

v velocity vector 

(x, y , z) moving coordinate system 
y surface-tension coefficient 

e small bookkeeping parameter 

f perturbation pressure 

X gas-to-liquid thermal diffusivity ratio 

p gas-to-liquid viscosity ratio 

p gas-to-liquid density ratio 

4> s perturbation in location of gas-liquid in- 

terface 

<P S location of gas-liquid interface 

co complex perturbation frequency 

Subscripts, Superscripts: 

i inner wave-number regime or integer vari- 

able 

f far outer wave-number regime 

l liquid 

g gas 

o outer wave-number regime 

* scaled quantity 
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Appendix. The Pulsating Stability Boundary 
for Small Liquid Prandlt Numbers 

For small liquid viscosities such that P = Pe and 
the same scalings introduced in the P ~ 0(1) case 
analyzed in the main body of the text, the effect of 
viscosity on the cellular boundary disappears at lead- 
ing order (P — > 0 in equations 39 and 40), while the 
effect of viscosity on the pulsating boundary is only 
significant in the far outer wave-number regime. In 
that case, the appropriate expansions analogous to 
those introduced for the P 0(1) case are given by 


ico 

~ £~ m (iCO {) 

+ icoiS U4 

+ i(O iy £ m 

+ • • *) and 


r - 



r ( - 1/2) “ 

-i<Ovp* 


<7 ' 


+ 0(tr‘), 

9<- V4> ” 

(ico { JP) m 

(60) 

b, 

= 2 
i 

+ 

1 

-7/4 




+ tt.-We- 
1 

M + . . ., 

i = 1, 2, 

8 

(61) 

b, 

- 2 

+ b\ v,l £ 

+ . . . 




i = 3 , 4 , 5 , 6 

b ( = &[-*' 4>e-5« + fcj-Vi + 7 (62) 

Substituting these expansions into equations 18-25 
and equating terms corresponding to like powers of 
e then gives, as previously, a sequence of equations 
for the recursive determination of the coefficients in 
the foregoing expansions. Similar to the calculation 
for P ~ 0(1) in the outer wave-number regime, we 
obtain in this case that 

(icoo) i 2 « 2 (k}/p*2)( A ; - A*), icO] =0, 

ico 2 = kj-[(A*/p*) 2 - 2 Pkf - 1] (63) 

where A* ~ (p*/2)(p*ykf - 1) is the inviscid cel- 
lular boundary in the far outer wave -number regime 
given by the last of equations 28. The first of equa- 
tions 63 thus recovers the cellular stability boundary, 
but because io) 0 is purely imaginary for A* < A*, 
stability in that region is determined by the real part 


of the next nontrivial coefficient in the above expan- 
sion for ico. Thus, setting ico 2 = 9 in the last of equa- 
tions 63, the pulsating stability boundary in the far 
outer wave-number regime is given by 

a; ~ -p* (i + 2 %) 1/2 (64) 

which, in the limit fy— > 0, matches with the leading- 
order pulsating boundary A* ~ ~ p* in the outer 
wave-number region, which is unaffected by viscos- 
ity to this order of approximation. Thus, equation 64, 
which in terms of k is given by A* ~ —p*( 1 + 
2p£k) l/2 y is valid for arbitrary wave numbers. Writing 
P in terms of its unsealed counterpart P, this ex- 
pression becomes identical to equation 59, which 
thus remains valid in the limit of small P . 

Acb lowledgments 

This work was supported by the U.S. Department of 
Energy under Contract DE- AC 04-94 AL85000 and by the 
NASA Microgravity Science Research Program under 
Contract C-32031-E. 


REFERENCES 

1. Landau, L. D., "On the Theoiy of Slow Combustion ” 
Acta Phtjsicochim. t/HSS 19:77-85 (1944); Zh. Eksp. 
Teor. Fiz. 14:240 (1944). 

2. Vosen, S. R., in Twenty-Second Symposium (Interna- 
tional) on Condnistion, The Combustion Institute, Pitts- 
burgh, 1989, pp. 1817-1825. 

3. Levich, V. G , Dokl. Akad. Nauk SSSR 109:975-978 
(1956). 

4. Armstrong, R. C. and Margolis, S. B., in Twenty- Second 
Symposium (International) on Combustion, The Com- 
bustion Institute, Pittsburgh, 1989, pp. 1807-1815. 

5. Armstrong, R. C. and Margolis, S. B., Combust. Flame 
77:123-138 (1989). 

6. Margolis, S. B., Combust. Flatne 113:406—423 (1998). 

7. Williams, F. A., Combustion Theory , 2nd ed., Benjamin/ 
Cummings, Menlo Park, CA, 1985. 

8. Bechtold, J. K. and Margolis, S. B., S/AM /. Appl. Math. 
51:1356-1379(1991). 


